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Introduction 


HULL  is  a  complex  ar  versatile  hydrodynamic  computer  code.  The 
original  version  was  dev*  -^ped  at  the  Air  Force  Weapons  Laboratory  in  1976 
[1].  It  has  the  capability  to  run  Eulerian  or  Lagrangian  modes  of 
calculations  and  has  a  wide  range  of  options  representing  material  selection 
and  various  physical  behavior.  One  of  the  drawbacks  of  HULL  is  its  explicit 
mode  of  calculation  with  a  severe  time-step  limitation.  The  time-step  of 
calculation,  restricted  by  the  Gourant,  Friedricks  and  Lewy  (GFL)  condition, 
often  hampers  its  efficiency  in  solving  large  multidimensional  problems.  Here 
size  and  complexity  of  the  code  adds  a  certain  amount  of  rigidity  in  upgrading 
HULL.  It  would  be  ideal  if  calculations  of  interest  to  the  Army  could  be 
speeded  up  substantially. 

A  prior  report  [2]  by  the  author  noted  certain  anomalies  along  the 
centerline  of  HULL  axisymmetric  calculations.  In  the  previous  study,  airblast 
velocity  profiles ,  pressure  plots ,  and  the  time  of  arrival  at  ground  zero , 
were  computed  from  an  airburst  detonation.  The  need  for  a  quicker  solution  in 
such  an  airblast  calculation  was  apparent.  The  present  effort  concentrates  on 
upgrading  HULL  by  adding  a  one -dimensional  spherical  airblast  code,  SAP 
(Spherical  Air  Puff). 

The  original  SAP  code  was  developed  at  the  Air  Force  Weapons  Laboratory 
in  the  1960's  [3].  However,  the  version  [4]  which  was  used  in  this  study  was 
developed  by  the  S-Cubed  Corporation  in  the  mid  1980‘s.  It  has  a  stand-alone 
capability  of  airblast  calculations.  SAP  can  also  prepare  station  plots  or 
feed  in  inputs  to  HULL  or  SHARC  calculations.  Unfortunately,  the  version  that 
was  available  to  the  author  lacked  several  subroutines.  The  need  for  complete 
documentation  resulted  in  redefining  several  system  parameters  and  options. 

The  program  was  added  to  the  HULL  library  in  UNICOS  and  VAX  modes.  However, 
thorough  tasting  of  SAP  against  measured  airblasts  still  needs  to  occur. 

This  report  presents  the  discussion  of  the  parameters  and  options  in  the 
newly  added  capability.  For  brevity,  only  the  important  details  will  be 
presented  here .  The  executable  versions  and  procedures  are  given  in  program 
listings  in  the  appendix . 


HULL  SYSTEM 


Before  the  discussion  of  SAP  begins,  it  is  important  to  recognize  and 
understand  the  HULL  framework.  The  HULL  system  consists  of  several  programs 
which  furnish  pre-  and  post-processing  needs  in  addition  to  the  HULL  program, 
which  performs  the  hydrodynamic  calculations .  The  programs  share  a  common 
library  of  utility  routines .  HULL  is  very  flexible  in  being  able  to  run  on 
several  computer  systems  and  can  be  upgraded  very  easily  because  of  its 
modular  structure.  HULL  is  prepared  by  a  proprietary  program  called  SAIL, 
developed  by  Shalimar  Research,  Inc  [5].  This  code  guides  the  flow  of 
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instructions  in  the  HULL  system  by  using  statements,  called  SAIL  directives, 
which  prepare  or  change  the  compilable  HULL.  This  is  important  to  recognize 
because  the  SAP  code  which  was  added  to  the  HULL  system  would  not  function 
itself  without  being  processed  by  SAIL. 

The  details  of  the  SAIL  program  cannot  be  published  because  of  the 
proprietary  guidelines.  It  operates  under  three  major  modes  - 
executive /maintenance ,  update/generate  and  copy/convert.  The  SAIL  maintenance 
and  executive  directives  are  its  normal  mode  of  operation.  However,  they 
cannot  be  processed  simultaneously  with  the  update/generate  mode.  Once  a 
program  is  transported  between  systems  using  the  copy/convert  mode,  it  can  be 
added  or  upgraded  for  another  installation  using  proper  SAIL  directives. 
Finally,  once  system  parameters  and  options  are  defined  and  available,  a 
program  can  be  executed  in  the  normal  mode. 

The  normal  mode  of  execution  is  prepared  for  HULL  with  the  pre¬ 
processing  programs  PLANK  and  KEEL.  The  preparation  of  the  compilable  HULL  is 
achieved  by  PLANK,  which  defines  the  options  block.  KEEL  utilizes  the 
material  library  and  generates  the  input/restart  tape  for  hydrodynamic 
calculations  of  HULL.  Finally,  there  are  post-processing  programs  such  as 
PULL  and  STATion  to  plot  the  results  of  hydrodynamic  calculations  and  analyze 
data. 


SAP  Code 


SAP  is  similar  to  the  program  HULL.  It  already  contains  post-processing 
programs  such  as  SAPPLT  and  STAPLT.  It  was  therefore  unnecessary  to  share 
those  programs  from  the  HULL  system.  Further,  SAP  does  not  need  the  HULL 
material  library  "MATLIB” ,  since  materials  are  furnished  in  several 
subprograms  within  SAP  which  can  be  compiled  upon  need.  Currently  SAP  has 
24  materials  in  its  library  (these  include  burned  and  unbumed  materials). 
They  are  -  AIR,  TNT,  BTNT,  PBX,  BPBX,  PENT,  BPENT,  ANFO,  BANFO,  METH,  BMETH, 
WATER,  HMX,  BHMX,  COMB,  BCOMB,  FOAM,  CFOUR,  BCFOUR,  PETN,  BPETN,  IRMT,  B1RMT 
and  PLUG.  The  material  names  beginning  with  the  letter  "B"  represent  burned 
explosives.  A  continuous  burn  routine,  based  upon  the  Chapman  Jouget  Theory, 
is  available  in  SAP. 

SAP  also  has  the  capability  of  using  six  different  ambient  models.  They 
are:  Tropical,  Temperate,  Arctic,  Exponential,  Constant  Energy  and  Hiroshima. 
SAP  calculations  can  incorporate  radiation  loss.  There  are  three  radiation 
loss  schemes  -  one  based  upon  SPUTTER  calculations,  one  using  emissivity  data 
by  Gilmer  [3],  and  a  radiation-diffusion  package  with  free  streaming  at  the 
edge  of  the  temperature  gradient.  SAP  uses  a  real  atmosphere  that  is  stable 
under  an  r'2  gravity  field,  and  an  equation  of  state  for  air  that  is  an 
empirical  fit  to  Hilsenrath's  data  [6], 

SAP  can  be  run  with  a  cartesian,  cylindrical  or  spherical  system  of 
coordinates.  It  has  two  differencing  schemes.  The  original  scheme  is 
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Lag  angian,  whereas  the  second  scheme,  which  is  also  utilized  in  the  HULL 
c,-ae,  and  can  be  run  as  Eulerian  or  Lagrangian,  as  set  by  an  option  SAP  can 
incorporate  artificial  viscosity  in  linear  or  quadratic  form.  It  even  has  a 
routine  for  an  artificial  compression  method.  The  SAP  plot  routines  can 
furnish  overpressures,  overdensity,  and  velocity  vectors  versus  radius.  SAP 
was  originally  developed  for  a  CDC  or  CRAY  machine  under  the  CTSS  mode . 
Therefore  to  run  under  CRAY-UNICOS  at  WES,  some  parameters  needed  to  be 
altered  and  some  new  ones  introduced  that  would  be  suitable  to  the  system. 


SAP  Options 


SAP  has  more  than  50  options  which  need  to  be  defined  to  prepare  a 
compilable  version.  If  one  prepared  a  PLANK  version  of  Sf.P,  that  program 
could  define  the  options  in  the  pre-processing  run.  The  main  SAP  run  can  use 
those  defaults,  or  modify  some  dynamic  options  during  the  executable 
processing.  Such  a  program  was  prepared  by  the  author  and  currently  resides 
in  the  SAP  directory  under  the  name  SPLANK.  However,  unlike  HULL,  fewer 
options  in  SAP  require  dynamic  alteration.  Therefore  one  could  skip  the  use 
of  SPLANK  and  directly  introduce  the  options  during  the  update/generate  mode 
of  SAIL. 

Some  options  also  have  dependency  on  other  options.  Therefore  care  must 
be  exercised  in  altering  them.  Appendix  A  shows  a  list  of  default  and 
redefined  options  as  listed  by  SAIL  during  pre-processing  of  SAP.  SAP  was 
added  to  the  HULL  system  as  the  procedure  RUNSAP.  Some  of  the  printed  options 
in  the  list  therefore  come  from  the  HULL  options  and  the  system  parameters. 
Appendix  B,  which  was  reproduced  from  [4],  gives  a  list  of  SAP  options  and 
overviews  of  the  subroutines  in  SAP.  One  must  keep  the  default  options  in 
mind  while  upgrading  the  SAP  system.  The  defaults  need  not  be  altered  at  all 
if  the  physical  problem  does  not  require  it.  In  the  present  CRAY  version  of 
SAP,  the  user  needs  to  cl  rose  the  following  options: 

1)  SAPNM  -  represents  the  number  of  materials  in  the  problem. 

2)  SAPMX  -  represents  the  maximum  grid  size. 

3)  SPGEOM  -  chooses  the  type  of  geometry  in  the  run. 

(Cartesian  -  1,  Cylindrical  -  2  and  Spherical  -3) 

4)  SPMTHD  -  chooses  between  differencing  techniques  in  the  solution  process. 

5)  NMAT  -  represents  the  number  of  material  boundaries  in  the  problem.  Note 
that  this  number  may  be  larger  than  SAPNM,  depending  on  the  arrangement  of 
materials . 

6)  RESTART  &  HYDRA  -  These  are  similar  options.  The  latter  should  be  chosen  to 
generate  a  restart  file,  whereas  the  former  should  be  chosen  when  a  restart 
file  already  exists. 

7)  BURN  -  1  must  be  used  if  explosives  are  present  in  the  problem.  In  those 
cases ,  the  explosives  can  have  both  burned  and  unburned  choices  (material  name 
preceded  by  "B"). 

8)  Other  important  options  include  RAD,  VISC,  ACM  and  FCT  to  specify  physical 
and  computational  choices,  whereas  SPUTIN,  HCUTIN  and  PLOTIT  give  pre-  and 
post-processing  choices. 
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SAP  INPUT 


After  the  pre-processing  SAIL  run  produces  a  compilable  SAP,  it  must  be 
compiled,  linked  and  executed  with  a  proper  SAP  input.  See  the  example  in 
Appendix  C.  SAP  inputs  have  generally  three  sections. 

The  first  section  starts  with  the  keyword  "sap” .  This  section  enters 
the  problem  number  and  starting  requests  in  terms  of  time  and  cycles.  The 
problem  number  and  starting  information  must  be  carefully  inputted  on  a 
restart  run,  the  failure  of  which  will  immediately  abort  the  run. 

The  second  section  starts  with  the  keyword  "input".  This  section  has 
the  stopping  request,  the  header  block  definition  and  definitions  of  the  "z- 
block"  parameters.  These  parameters  are  upgraded  at  each  run  and  therefore 
some  of  these  can  be  considered  as  dynamic  options.  A  complete  list  of  the  z- 
bloek  parameters  are  available  in  appendix  D. 

The  last  section  starts  with  the  keyword  "region".  This  is  similar  to 
the  keyword  "mesh"  in  KEEL  runs.  It  has  six  different  variables  separated  by 
the  word  CELLS.  These  variables  are  -  material  name,  dx,  energy,  rho,  speed, 
and  bound.  Material  names  have  to  be  used  exactly  as  they  were  introduced  in 
the  SAIL  pre-processing.  DX  represents  the  width  of  cells  in  this  part  of  the 
mesh . 


Energy  and  rho  represent  the  energy  and  mass  densities.  All  units  are 
expressed  in  the  CCS  (Centimeters,  Grams,  Seconds)  unit  system.  Speed 
represents  the  speed  of  the  material  in  this  part  of  the  mesh.  Again,  this 
quantity  is  expected  to  be  a  variable.  The  first  cell  must  be  assigned  some 
initial  velocity.  A  zero  velocity  input  for  the  first  cell  will  result  in 
aborting  the  program.  Bound  represents  the  maximum  number  of  cells  for  which 
the  previous  five  variables  hold.  The  default  on  bound  is  IMAX,  the  maximum 
number  of  cells  in  the  problem.  If  there  are  several  materials  present  in  the 
problem,  the  bound  values  should  be  inputted  carefully. 

All  of  the  above  six  variables  may  not  be  specified  in  a  given  problem. 
For  example,  the  sample  input  in  Appendix  C  has  three  materials  -  SPENT,  PENT 
and  AIR.  The  problem  stacks  up  BPENT  in  the  first  cell,  PENT  in  Cell  2 
through  77  and  AIR  for  Cells  78  through  100.  The  mesh  size  is  a  variable.  As 
may  be  noticed  from  this  example,  only  4  variables  (out  of  6)  were  defined  in 
the  first  three  cells,  3  were  defined  between  Cells  4  through  77  and  only  two 
were  defined  for  the  rest  of  the  mesh.  Further  descriptions  of  input 
variables  are  presented  in  Appendix  D. 

The  final  word  in  the  sap  input  is  "end".  Once  these  inputs  are 
correctly  specified,  SAP  outputs  a  "sap. out"  and  a  "tape  4"  (if  option  hydra 
is  chosen) . 
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Results 


Appendix  E  shows  the  sequence  of  the  batch  run  "jobsap.com".  This  batch 
job  prepares  the  compilable  version  from  the  card  input  [4],  compiles,  links, 
and  executes  the  program  with  the  SAP  input  shown  in  Appendix  C.  The  output 
"sap. out"  is  presented  in  Appendix  F.  It  can  be  seen  from  the  output  that  the 
program  generates  the  2 -block  and  the  mesh  configuration  correctly. 

This  particular  run  ha^  a  problem  in  the  subroutine  "BURN",  as  reflected 
in  the  output.  This  code,  like  any  other  newly  developed  capability,  has  to 
be  tested  in  various  physical  problems  before  confident  statements  can  be  made 
regarding  its  advantages  and  drawbacks .  The  error  source  is  probably  among 
input  or  material  properties  data  or  methodology  of  calculations.  Further 
testings  are  currently  in  progress. 


Conclusions  and  Recommendations 


The  SAP  code  was  added  to  the  HULL  library  of  programs  at  WES  This 
code  will  provide  an  input  for  the  HULL  calculations,  or  it  can  produce  one- 
dimensional  solutions  for  airblast  problems  independently.  The  code  remains 
to  be  tested  in  various  physical  problems. 

One  very  useful  application  of  the  code  will  be  to  use  it  in  the 
airblast  calculations  presented  in  [2].  It  would  be  interesting  to  know 
whether  SAP  outputs  can  be  inputted  into  HULL  calculations  with  significant 
improvement  in  the  overall  run  time  of  the  HULL  calculations .  Another 
interest  would  be  to  see  if  the  HULL  input  thus  generated  by  the  "firein" 
option  in  KEEL  would  produce  a  better  arrival  time  of  the  shock  at  ground  zero 
in  the  HULL  run  than  was  obtained  previously. 

SAP  can  also  be  used  in  the  reverse  mode ,  where  HULL  outputs  are  used  as 
input  to  SAP  by  the  program  HULLCTJT,  which  takes  a  slice  of  the 
multidimensional  output.  Finally,  the  newly-added  capability  of  one- 
dimensional  calculations  can  produce  approximate  and  quick  answers  to 
significantly  large  problems  as  a  stand-alone  capability,  where  an  early 
prediction  can  influence  later  decision  making. 
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Appendix  K 


1  sail  library  for  version  124  of  system  hull 

generated  3aug92 


this  library  contains  140358  lines 
generating  source  from  programs  on  3aug92 
program  origin  program  origin 

program  origin 


*sap  210000 

programs  which  are  marked  with  an  *  will  generate  source  in 
this  run 

1  options  defined  at  the  beginning  of  this  run 


0  inst  = 

0  runsap  = 

0  matlib  = 

0  proc-msc  = 

0  pull  = 

0  keel  = 

0  staplt  = 

0  library  = 


16  mach  = 
1  tekcolor  = 
0  plank  = 
0  hull  = 
0  station  = 
0  lrez  = 
0  hullcut  = 
0  control  = 


loptions  defined/redefined  during 
cray  =  1 

opsys  =  2 

vendor  =  1 

pf  =  1 

ob j lib  =  0 

tape lib  =  3 

plotpkg  =  0 

remark  =  0 

mf  =0 

cdc  =  0 

ibm  -  0 

cray  =  1 

vax  =  0 

sel  =  0 

vsos  =  0 

colour  =  0 

word  =  64 

byte  =  8 

date  =  2 

buffer  =  1 

sensesw  =  1 

nw  *=  1 

CW  ®=  8 

cardl  *=  10 


1  sail  = 

0  eos  = 

0  common  = 

0  plotters  = 

0  plotvg  - 

0  sap  = 

0  sapplt  = 

0 

executive  processing 


0 

0 

0 

0 

0 

2 

0 
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cardo  - 

double  = 

exbl  “ 

debugl  = 

debugs  = 

mach  = 

cdc  = 

univac  - 

vax  = 

cel  - 

ridge  = 

spsvth  = 

sapmx  = 

spgeom  - 

spmthd  = 

sapnm  = 

single  = 

airopt  = 

atmos  = 

compare  = 

nstat  = 

nmat  = 

matt  -- 

air  = 

tnt  = 

btnt  = 

pbx  = 

bpbx  = 

pent  = 

bpent  = 

anfo  " 

banfo  = 

meth  = 

bxneth  - 

water  = 

hmx  = 

bhmx  * 

comb  ~ 

bcoiab  ~ 

foam  ~ 

cfour  «* 

bcfour  *= 

petn  = 

bpetn  = 

irmt  = 

tairmt  = 

plug  = 

rezone  *= 

method  ~- 

lbnd  ~ 

rbnd  ~ 

imax  = 


8 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

100 

3 

1 

3 

0 

1 

2 

0 

0 

3 

24 

3 

0 

0 

0 

0 

2 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

1 

0 

0 

100 


OHO 


restart 

ss 

0 

ktstd 

— 

0 

weight 

= 

1 

ngroup 

ass 

0 

sputin 

*= 

0 

hcutin 

ss 

0 

hydra 

St 

1 

acm 

0 

burn 

SS 

1 

fct 

= 

0 

rad 

ss 

0 

run 

ss 

3 

plotit 

== 

0 

geom 

ss 

3 

vise 

ss 

1 

euler 

SS 

0 

jwl 

ss 

1 

cray 

ss 

1 

unicos 

ss 

1 

eos 

= 

1 

sum 

ss 

1 

route 

ss 

0 

sys 

= 

2 

nchpwd 

ss 

s 

nodpwd 

— 

22 

nbitpc 

ss 

8 

nchar 

= 

16 

lev 

3= 

1 

vector 

SS 

0 

asciim 

s 

1 

disspl 

* 

0 

scaler 

= 

1 

nvarpr 

SS 

500 

nh 

SS 

5 

nbitpw 

= 

64 

nshift 

SS 

32 

nshif 1 

ss 

56 

nchmtn 

= 

2 

nwpcv 

ss 

2 

nzblk 

ss 

100 

indxz 

ss 

200 

nzblk2 

ss 

200 

idim 

ass 

300 

ijump 

ss 

101 

ijptwo 

= 

201 

2856  cards  generated 

end  of  normal  sail  run 

system  hull  version 


124 
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Appendix  B 

*p 

options 

20feb87 

20feb87 

as 

option  < 

default 

definition 

20feb87 

20feb87 

ss 

atmos 

2 

atmosphere  type 

20feb87 

ss: 

20feb87 

=1  tropical  atmosphere 

2E 

20feb87 

=2  temperate  atmosphere 

20feb87 

=3  arctic  atmosphere 

zr 

20feb87 

=4  exponential  atmosphere 

3= 

20feb87 

=5  constant  energy  and  pressure 

* 

20feb87 

=6  hiroshima  atmosphere 

= 

air opt 

1 

air  type,  this  option  is  for  use  in 

20feb87 

20feb87 

conjunction  with  vector=l,  cray=l 

= 

20feb87 

=1  doan-nickel  equation  of  state 

= 

20feb87 

=2  gamma  law  equation  of  state 

20feb87 

=3  simple  gamma  law  equation  of  state 

= 

20feb87 

=10  four  point  interpolation 

as 

vector 

0 

vectorized  eularian  differencing,  to  be 

used 

20feb87 

= 

20feb87 

only  with  cray=l 

s 

20feb87 

=0  no  vectorization 

as 

20feb87 

=1  use  vector  routines 

as 

Cray 

0 

sap  run  to  be  run  on  cray-1  computer 

20feb87 

= 

20feb87 

=0  not  a  cray  run 

20feb87 

=1  cray  run 

1 

i _ 

10 

r/  • 


=  cdc 

20fcb87 

1 

sap  run  to  be  run  of  cdc  computer 

20feb87 

~0  not  a  cdc  sap  run 

20feb87 

=1  cdc  sap  run 

=  restart 

20feb87 

0 

this  option  to  restart  a  previous  run 

20feb87 

=0  not  a  restart  run 

20feb87 

=1  problem  is  a  restart  run 

=  euler 

20feb87 

0 

reference  frame  option 

20feb87 

=0  lagrangian 

20fefc87 

=1  eulerian  {  use  only  with  method  2 

20feb87 

single-material  runs.) 

-  geom 

20feb87 

3 

geometry  option 

20feb87 

=1  cartesian 

=2  cylindrical 

20fab87 

=3  spherical 

=  Ktstd 

0  1  kiloton  standard  left  bdry  input 

condition  20feb87 

20feb87 

=0  not  included 

20feb87 

=1  included,  method  2  should  be  used. 

20feb87 

subroutines  as  contained  on  sail  file 

20feb87 

"dytaid"  under  the  name  Iktstd  should  be 

20f..bS7 

loaded  with  sap  if  this  option  is  used. 

this  20f«b87 

options  lbnd  and/or  rbnd  must  be  set  if 

*S 

is  to  be  used. 

20ftb87 

=  Xbnd 

20feb87 

0 

left  boundary  condition  option 

20feb87 

=0  no  left  boundary  input 

■  imposed.  2Gfeb87 

=1  there  is  a  left  boundary  condition 

called  2Qfeb87 

if  ktstd=0„  a  left  boundary  subroutine 

f 

11 

t 


4 


sb  " Ibound (time , distance , u ,  rho ,  i )  M 

(method  1)  20feb87 

=  " Ibound (time, u, rho, i) " 

(method  2)  20feb87 

=  which  specifies  speed,  density  and  energy 

as  2Qfeb87 

-  a  function  of  time  (and  distance,  if 

method  1)  2Gfeb87 

"  must  be  loaded. 

20feb07 

*  also,  the  parameter  "lbref"  set  by  input 

must  20feb87 

=  be  zero,  indicating  a  transmissive  left 

2G£@b87 

-  boundary 
20feb87 

=  rbnd  0  right  boundary  condition  option  (method  2 

only)  20feb87 

-*  -0  right  boundary  condition  is  defined  bv 

20feb87 

"  the  option  ' rezone' 

20feb87 

“  =1  the  right  boundary  is  reflective,  the 

20feb87 

=  sap  variable  ' shkcel '  will  never  be 

permitted  20feb87 

~  to  have  values  greater  than  imax  -  2,  and 

the  20feb87 

=  calculation  will  not  terminate  when  shkcel 

20feb87 

~  reaches  a  value  of  imax  -  2. 

2Gfeb87 

~  =2  there  is  a  right  boundary  condition 

imposed  20feb87 

=  if  ktstd-O,  a  right  boundary  subroutine 

called  20feb87 


*  "r bound ( time , u , rho , i ) " 

(method  2)  20feb87 

-  which  specifies  speed,  density  and  energy 


as 

20feb87 

ss 

20feb87 

a  function  of  time  must  be  loaded. 

=5 

method 

2 

differencing  scheme  option 

20feh87 

20feb87 

=1  old  sap  differencing 

20feb87 

—2  hull  differencing 

ss' 

weight 

1 

weight  option  for  method  2 

20feb87 

ss 

20feb«7 

=1  conventional  hull  boundary  algorithm 

12 


(  bv=(vl+v2)/2  ) 


20feb87 


20feb87 


™2  mass  weighting  boundary  algorithm 


20feb87 

=  nstat 

20feb87 
-  nmat 

20febS7 

20feb87 

~  rezone 

2Gfeb87 

20feb87 

while  20feb87 


energy)  20feb87 
while  20feb87 
20feb87 

s z 

20feb87 

=  nvarpr 

20feb87 

20feb87 

20feb87 

20feb87 
=  imax 

20feb87 
=  rad 

20feb87 


20feb87 


sputter 


20feb87 


20febS7 


20ieb87 


0  number  of  stations  in  problem 

1  number  of  material  boundaries  in 

problem  at  any  given  time  (see  below) . 

1  rezone  option 

=0  no  rezone  (stop  at  right  edge  of  mesh) 
=1  rezone  (expand  front  fourth  of  mesh 
conserving  mass,  momentum  and  total 
=2  rezone  (expand  front  fourth  of  mesh 
conserving  mass,  momentum  and 
mass*pressure) 

400  maximum  number  of  hydro  variables  per 

tape4  hydro  record 

=imax*5  for  method  1 

= imax* 4  for  method  2 

100  number  of  cells  in  problem 

0  radiation  option 

=0  no  radiation  loss 

=1  prescribed  power  out  based  upon 

r 

calculations 

=2  emissivity  based  upon  data  by  gilmer 


free-streaming  at  20feb87 


20feb87 


20feb87 


=3  radiation  diffusion  with 
the  edge  of  the  temperature  gradient. 


=4  -hfb  must  be  set  if  rad=4-  radiative 
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20feb87 

20feb87 

20feb87 

=  ngroup  1 

(used  20feb87 

20feb87 

*  sputin  0 

20feb87 

tape  20feb87 

tape  20feb87 

20feb87 

20feb87 

20feb87 

=  hcutin  0 

20feb87 

20feb87 

20feb87 

tape4  20feb87 

contains  20feb87 
and  20feb87 
if  20feb87 

can  20feb87 

20feb87 
20feb87 

=  plotit  0 

20feb87 

20feb87 

over-  20feb87 
a  20feb87 

20feb87 


heating  using  rate  dependant  thermal 
layer  modal .  currently  used  only  when 
angle=90. 

number  of  frequency  groups  in  problem 
only  when  rad=3 ) . 
sputter  input  for  mesh 
=0  input  is  by  means  other  than  a  sputter 
=1  input  to  fill  mesh  is  from  a  sputter 
this  option  can  be  used  at  any  angle  with 
any  geometry. 

input  is  from  a  tape4  generated  by  program 
hullcut. 

=0  input  is  by  means  other  than  hullcut. 
=1  input  is  for  fill  the  mesh  form  a 
generated  by  hullcut.  the  tape4 
a  z-block,  a  'x'  coordinate  record, 
the  hydro  data,  the  tape  is  used  as 
it  were  a  restart  tape,  the  option 
used  at  any  angle  with  any  geometry. 

plot  option 

=0  no  plot  file  is  generated 
=1  tapel  metafile  is  generated  to  plot 
pressure,  overdensity  and  velocity  as 
function  of  time. 
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3= 

hydra 

1 

20feb87 

20feb87 

20feb87 

2= 

20£ebS7 

SZ 

vise 

1 

20feb87 

20f eb87 

20feb87 

20feb87 

= 

a  cm 

0 

20febS7 

20feb87 

mass. 

20feb87 

compression  20feb87 


20feb87 

=  fct  0 

20feb87 

20feb87 

=  burn  0 

2Cfeb87 

sphere)  20feb87 


20feb87 

nm  1 

20feb87 

matt  24 

20feb87 

air  1 

20feb87 

tnt  0 

20fab87 

btnt  ~  0 

20feb87 

pbx  0 

20feb87 

bpbx  0 

20feb87 

pent  0 

20feb87 


restart  option 

=0  no  tape  4  is  used,  i.e,,  no  restart 
capability 

=1  a  tape4  (restart  tape)  is  written 

artificial  viscosity  option 

=0  no  artificial  vise  terms  added 

=1  constant  term  added 

=2  variable  term  added 

artificial  compression  method  option 

=0  no  artificial  compression  used. 

=1  artificial  compression  applied  to 
momentum  and  total  energy  in 
regions. 

flux  corrected  transport  scheme 


high  explosive  material  option 
=0  no  material  is  present  (isothermal 
=1  explosive  material  is  present 
number  of  materials  in  the  problem 
number  of  materials  available  in  sap 
“positive  value  if  present  in  problem 
“positive  value  if  present  in  problem 
^positive  value  if  present  in  problem 
“positive  value  if  present  in  problem 
“positive  value  if  present  in  problem 
“positive  value  if  present  in  problem 
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= 

bpent 

20feb87 

0 

=positive 

value 

if 

present 

in 

problem 

meth 

20feb87 

0 

“positive 

value 

if 

present 

in 

problem 

ss 

bmeth 
20f eb87 

0 

“positive 

value 

if 

present 

in 

problem 

anfo 

20feb87 

0 

=positive 

value 

if 

present 

in 

problem 

ss 

banfo 

20feb87 

0 

=positive 

value 

if 

present 

in 

problem 

=  water 

20feb87 

instructions) 

0 

20feb87 

=positive  value  if  present 

(refer  to  subrout i 

in 

ne 

problem 

waters  : 

= 

hmx 

20feb87 

0 

=positive 

value 

if 

present 

in 

problem 

sr 

bhmx 

20feb87 

0 

“positive 

value 

if 

present 

in 

problem 

ss 

comb 

20feb87 

0 

“positive 

value 

if 

present 

in 

problem 

ss 

bcomb 

20feb87 

0 

“positive 

value 

if 

present 

in 

problem 

ss 

foam 

20feb87 

0 

“positive 

value 

if 

present 

in 

problem 

= 

cfour 
20f eb87 

0 

“positive 

value 

if 

present 

in 

problem 

= 

bcfour 

20feb87 

0 

“positive 

value 

if 

present 

in 

problem 

ss 

petn 

20feb87 

0 

“positive 

value 

if 

present 

in 

problem 

bpetn 

20feb87 

0 

“positive 

value 

if 

present 

in 

problem 

= 

irmt 

20feb87 

0 

“positive 

value 

if 

present 

in 

problem 

birmt 

20feb87 

0 

“positive 

value 

if 

present 

in 

problem 

= 

plug 

20feb87 

20feb87 

0 

“positive 

value 

if 

present 

in 

problem 

the  integers  from  1  to  na  should  be  assigned  to 
20feb87 

materials  in  the  problem,  all  other  materials  should 
20feb87 

be  set  to  zero. 

20feh87 

2Gfeb87 

a )  a ir ( isotherma 1  sphere ) 

20feb87 

b)  bpent,air 
20feb87 


16 


c)  bpent, pent, air 

20feb87 

d)  water, bpent, pent, water 

20fab87 

e)  bpent, pent, tnt, pent, air 

20f eb87 

20feb87 

should  have  the  following  options  set  as  shown- 


20feb87 

20feb87 

20feb87 

20feb87 

nm 

nmat 

matt 

burn 

(material  codes) 

a) 

20feb87 

1 

1 

12 

0 

air=l 

b) 

20feb87 

2 

2 

12 

0 

air =2 , bpent=i 

c) 

20feb87 

3 

3 

12 

1 

bpent=l , pent=2 , air=3 

d) 

20feb87 

3 

4 

12 

1 

water=l , bpent=2 , pent=3 

e) 

20feb87 

20feb87 

4 

5 

12 

1 

bpent=l , pent=2 , tnt=3 , air=4 

except  for  a)  (obviously)  the  ordering  of  the 
20feb87 

material  code  numbers  is  arbitrary  (the  fifth  column 
20feb87 

of  the  above  example  is  not  unique) .  note  that  the 
20feb87 

last  zone  in  the  mesh  must  be  included  as  a  material 
20feb87 

boundary . 

20feb87 

20feb87 

sap  consists  of  a  fair  number  of  subroutines  linked  by 
20feb87 

blank  and  labeled  common,  in  addition  to  the  summary 
20feb37 

given  below,  most  of  the  subroutines  contain  further 
20feb87 

explanation  of  their  functions,  for  a  detailed 
20feb87 

description  of  card  input  to  sap,  see  subroutine  in. 
20feb87 

the  section  labeled  "constant"  in  sap  outlines  the 
20feb87 


procedure  for  adding  new  materials  to  sap. 

20feb87 

20febS7 

20feb87 

initializing  routines  -  these  are  called  when  setting  up 
20feb87 

or  restarting  a  sap  problem,  they  are  not  contained 
20feb87 

within  the  main  sap  loop,  an  asterisk  before  the  routine 
20feb87 

indicates  that  the  routine  is  included  only  if  the  option 
20feb87 

parentheses  after  the  description  is  true. 

20feb87 


20feb87 

in~in2  card  input  to  sap  (see  in  for  explanation  of 
20feb87 

input  format) 

20feb87 

♦restart  tape  input  (tape4 , tape9)  restart 
20feb87 

if  applicable  (hydra) 

20feb87 

cell  sets  up  mesh  using  input 

20feb87 

set  defines  zblk  default  values  before  input  is  read 

20feb87 

rstzbk  places  input  that  is  read  into  the  zblk  (from 
20feb87 

tape  or  card  input)  into  appropriate  problem 

20feb87 

parameters 

20feb87 

out  prints  the  initial  zblk  and  the  problem 

20feb87 

conf igurat ion . 

20feb87 

20feb87 

main  loop  routines 
20feb87 


20feb87 

dtime 

20feb87 

hydro 

20feb87 

20feb87 


selects  time  step  for  next  cycle 
calculates  hydro  variable  values  throughout 
active  mesh  for  the  current  cycle. 
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lbl  defines  left  boundary  values  for  method  2 

20f efoS 7 

lb3  defines  left  bdry  fluxes  for  eulerian  method  2 

20feb87 

*burn  advances  detonation  front,  converts  h.e. 

20febS7 

materials  and  adds  chemical  energy  released  by 

20£eb87 

burning  to  the  burning  cell,  (burn) 

20£ebS7 

*radlos/fblossl  energy  loss  by  radiation,  routine 
20feb87 

included  depends  on  sail  option  value,  (rad) 


20feb87 

*rezone 

20feb87 

*stuff 

20feb87 

20feb87 


redefines  mesh  coordinates,  (rezone) 
redistributes  hydro  quantities  in  the 
mesh  set  up  by  subroutine  rezone,  (rezone) 


20febS7 


=  utility 

20feb87 

20febS7 

=  *decide 

20feb87 

gtl)  20feb87 
»  *const 

20feb87 

20feb87 
=  *eos 

20feb87 

20feb87 
=  *visc2 

20feb87 

(visc2)  20feb87 

=  *q 

to  20feb87 
20feb87 

=  *acmsub 

terms  20feb87 

20feb87 

20feb87 


routines 

selects  material  content  of  a  cell  according  to 
the  cell-s  position  wrt  the  boundary  array,  (nm 
selects  detonation  energy  and  initial  velocity 
of  detonation  for  h.e.  materials,  (burn) 
calls  the  appropriate  equation  of  state 
routine  for  the  cell  in  question,  (nm  gt  1) 
determines  the  magnitude  of  the  linear  viscosity 
coefficient  as  a  function  of  maximum  density, 
determines  the  artificial  viscosity  to  be  added 
pressure  term  in  the  hydro,  (vise  lt2) 
modifies  the  mass,  momentum  and  total  energy 
in  regions  of  compression  tc  preserve 
discontinuities.  (acm) 
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kill 
2 Of ebS7 

20febS7 

20feb87 

*rated 

20feb87 

20feb87 

20feb87 

*dairt 

20£eb87 

20f eb87 

20feb87 

20feb87 

output 

20feb87 

20feb87 
edit 
20feb87 
*wtape4 
20f eb87 
setz 
20feb87 
*draw 
20feb87 

20feb87 
*wtape9 
20f eb87 

20feb87 

20£eb87 

20feb87 

20feb87 


stops  the  program  on  a  model  condition  for  a 
small  number  of  recognized  foul  ups  in  the 
other  routines. 

thermal  heating  routine,  returns  ground 
level  temperatures  as  a  function  of  incident 
thermal  flux  for  a  given  time. (rad4) 
the  derivative  of  subroutine  airt.  returns 
the  derivative  of  temperature  wrt  energy, 
used  in  thermal  heating  calculations. (rad4) 

routines  -  most  called  at  std  times 

writes  to  file  output, 
hydro  tape4  output  (hydra) 
updates  the  zblk  array. 

plots  of  velocity,  overdensity  and  overpressure 
vs  radius  are  produced,  (plotit) 
station  tapes  output  (called  from  subroutine 
hydro  when  a  station  is  encountered)  (nstat) 
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Appendix  C 


sap  prcb^lll . 0  t-1 . e-4 

cycle-0,0  input  cstop=2100  detcel=1.0  rdw=ll . 0 
angle-0.0  head 

sap  run  pentolite  100  cells  29  july  92 
hob-0 . 0  >:0=0 . 0 
region 

bpent=l  speed=.01  dx=11.0  bound=l  cells 
speed=.01  pent=2  dx=10.  bound=3  cells 
pent=2  dx=10.0  bound=75  cells 
air=3  dx=30.  bound=77  cells 
air=3  dx=42 .  cells 
end 
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Appendix  D 


subroutine  in 
20febB7 

c 

20feb87 

0**************************************************************** 
*******  20feb87 

c* 

*  20febS7 

c*  this  subroutine  reads  the  card  input  to  sap  and  places  the 

*  20feb87 

c*  values  in  the  proper  parameters, 

*  20feb87 
c* 

*  20feb87 

^************ ********************************************** ****** 
*******  20feb87 
c 

20feb87 

$  parameter  (nmat  =  _nmat_  ) 

Ilf ebS7 

equivalence  (what,  iws) 

20£eb87 

20f ebS7 

~  input  to  sap  is  divided  into  several  sections  separated 

20f eb87 

-  by  keywords  given  below,  all  input  is  in  free  format. 

20feb87 

20feb87 

--  the  first  keyword  is  "sap”  and  must  in  all  cases  be  the 

20feb87 

=  first  word  of  data,  three  parameters  follow 
20feb87 

t  «=  (time  for  start  of  this  run.  if  not  present,  default 

20£eb87 

-  is  to  a  large  value,  indicating  a  restart  at  the  end  of 
20£eb87 

*  the  last  run.  the  time  may  not  be  . le.  zero,  if  input  as 

2Gfeb87 

=  such  the  program  will  stop) 

20feb87 

~  prob  *=  (problem  number) 

20feb87 

--  cycle  *=  (cycle  number  for  start  of  this  run.  for  a  new 

20feb87 

as  problem,  cycle*=0.0.  default  is  to  a  large 
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value, indicating  20feh87 

=  a  restart  at  the  end  of  the  last  run.  any  other  positive 

20£ebS7 

=  number  will  result  in  a  restart  at  that  cycle  number) 

20£eb87 

-  for  a  new  problem,  time  and  cycle  must  both  be  specified. 
20feb87 

=  for  a  restart  neither  need  be  input,  failure  to  specify 

20feb87 

=  the  correct  problem  number  on  a  restart  will  cause  sap  to 
20feb87 
=  abort . 

20£eb87 

20feb87 

=  for  the  previous  two  sections,  the  data  input  may  appear 

more  20feb87 

-  than  once.  the  latest  value  will  be  the  one  used  by  the 
program.  20feb87 

-  this  is  convenient  when  running  two  or  more  problems  which 
differ  20feb87 

-  only  in  the  values  of  the  input  parameters,  input  cards  may 

be  20feb87 

~  reordered,  rather  than  redone. 

20feb87 

20feb87 

=  the  next  section  is  delimited  from  the  first  by  the  word 

20feb87 

*  "input",  data  contained  in  this  section  pertain  to  zblock 

20feb87 

=  parameter  values,  header  card  information,  and  the  station 

20feb87 

-  coordinates  if  applicable,  with  the  exception  of  the  header 
20feb87 

=  card,  any  values  in  this  section  may  be  defaulted  to 

internally  20feb87 

=  set  values  listed  below,  data  in  this  section  may  appear  in 
20feb87 

-  any  order  with  two  exceptions,  the  header  card  must  be 
20feb87 

»  on  a  separate  card,  preceded  by  the  word  "head",  and  the 

20£eb87 

-  station  coordinates  or  data  for  generation  of  stations  must 

20feb87 

-  be  preceded  by  the  word  "stations"  (see  below) . 

20feb87 

20£eb87 

=  word  default  definition 

2C£eb87 


23 


20f eb87 

=  angle  0.0 

up  20feb87 
=  estop  l.el7 

20feb87 

-  cO  1.80 

20£ebS7 

=  cl  .5 

20feb87 

=  detcel  0 . 

20feb87 

regions.  20feb87 
=  dstop  l . e9 

20feb87 

=  dt  1.  e-9 

20feb87 

XS 

20feb87 

“  dxmul  1 . 0 

20feb87 

20feb87 

=  hob  0 . 0 

20feb87 

=  lbref  l 

20feb87 

20feb87 

20feb87 

20feb87 

20feb87 

=  prdw  0 . 0 
20feb87 

20feb87 

*  ptstop  100. 

20feb87 

=  qmin  . 1 

20£eb87 

=  rdw  0 . 0 

20feb87 

3 S 

sphere)  20feb87 

=  stabf  0 . 8 

20febS7 

20feb87 

=  umin  1 . e-4 


angle  of  run  (degrees)  +90  degrees  -  straight 
cycle  stop 

quadratic  artificial  viscosity  coefficient 
linear  artificial  viscosity  coefficient 
number  of  cell  containing  detonation  wave 
this  cell  separates  burned  from  unburned 
distance  of  max  pressure  stop  (cm) 
time  increment  (input  value  used  only  for 
first  cycle) 

increment  ratio  in  exterior  cell  size  (used 
in  rezone) 

height  of  burst  (cm) 
left  boundary  reflective  indicator 
(  a  value  of  0  -  indicating  transmissive 
left  bdry  -  is  sensible  only  with  method  2 
and  cartesian  geometry,  it  must  be  used  if  a 
left  boundary  input  condition  is  used.) 
radius  of  detonation  wave  last  cycle  (=0.0  for 
new  problem  -  set  by  tape4  read  on  restart) 
problem  time  stop  (sec) 

minimum  artificial  viscosity  (noise  level) 
radius  detonation  wave  this  cycle  (=0.0  for 
for  a  problem  with  no  burn  -  isothermal 
stability  factor  (safety  factor  in  courant 
condition  time  step) 
minimum  hydro  velocity  (noise  limit) 
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20febS7 

=  XO  0.0 

non-zei'o  20feb87 

20feb87 

=  yield  0.0 

20£abf»7 

=  hfb  hob 

20feb87 

=  grndr  0 . 0 

(rad4)  20feb87 


distance  of  mesh  left  edge  from  origin.  a 
value  makes  sense  only  if  geom=l.  (cm) 
yield  of  device  in  kt 

height  of  fireball  above  ground  (cm)  -  (rad4) 
distance  from  mesh  to  ground  zero  (cm)  - 


20feb87 

=  any  of  the  above  may  be  rstzbk  by  specifying  the  keyword 

given  20feb87 

=  followed  by  the  desired  value,  any  such  changes  must  appear 
20feb87 


in  the  section  labeled  ,(l  input". 

20feb87 

some  notes  on  the  keywords  -estop-,  -dstop-,  and  -pstop-,  if 
20feb87 

any  of  the  above  are  found,  then  the  ones  found  will  be  the 
20feb87 

only  conditions  which  will  stop  sap  during  a  run.  if  the 
20feb87 


=  user  omits  all  of  the  stop  conditions,  then  all  three 

default  20feb87 


values  will  be  used  for  ending  sap. 
20feb87 


20feb87 

=  the  sputin  section  should  follow  the  input  section,  and 

20feb87 

=  it  can  be  before  or  after  the  header  card  sequence.  the 

default  20feb87 

=  values  found  here  are  different  once  the  key  word  -sputin- 

is  20feb87 

=  encountered,  following  the  afore  mentioned  keyword  are  the 

20feb87 

=  following  parameters 

20feb87 


20feb87 

word  default  definition 

20feb87 


20feb87 
yield  0.0 
20feb87 


scale  n/a 
20feb87 


yield  of  device  in  kt. 

if  this  word  is  encountered,  then  the  data  is 
scaled  to  user  input  -  i.e.  yield  and  hob. 
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20f eb87 


ssz 

hot  0 . 0 

20feb87 

height  of  burst  in  cm. 

rs 

fb  none 

mandatory  name  of  sputter  tape  to  be  used. 

the 

20feb87 

XT 

actual  tape  itself  should  be  called  tape!  in 

the 

20feb87 

«»** 

20feb87 

users  local  file  space. 

*= 

dx  100. 

cell  width  after  the  sputter  tape  input  has 

been 

20feb87 

= 

put  into  the  mesh. 

20feb87 

20feb87 

*=  note  that  the  parameters  yield,  hob,  and  d>:  are  found  in 

other  20feb87 

-  input  sections,  these  values  must  appear  in  the  -sputin- 
2Gfeb87 

=  section  or  the  values  are  not  guaranteed  to  be  of  value. 

20feb87 

*=  it  is  also  important  to  note  that  the  -region- 

20feb87 

~  section  is  ignored  when  sputter  tapes  are  used  as  input,  if 

the  20feb87 

=  angle  in  the  input  section  is  >  0  and  if  the  -scale- 

parameter  is  20feb87 

=  set,  then  it  will  scale  the  data  to  a  mesh  cut  at  that 

angle.  20feb87 

20feb87 

-  stations  may  be  specified  if  desired  in  one  of  two  ways 
20feb87 

20feb87 

=  stations  explicit  (nstat  values  listed  in  increasing 

20f ebd  7 

=  crder,  to  be  used  as  coordinates  of  stations) 

20feb87 

=  -  or  - 

20feb87 

=  stations  xsO  =  (coordinate  of  first  station) 

20feb87 

=  dx  «  (distance  between  adjacent  stations) 

20feb87 

=  either  parameter  may  be  specified  first,  and  both  must  be 

present  20feb87 

*=  if  sail  option  "nstat”  is  set  to  0,  of  course  no  keyword  or 

20feb87 

=  data  pertaining  to  stations  should  be  present  in  input. 
20feb87 
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20f eb87 

the  final  section  of  data  is  begun  with  the  word  “region”. 
20feb87 

this  section  need  not  be  included  for  a  new  problem  if  option 
20£eb87 

•*sputinM  is  set. 

20feb87 

this  section  is  used  to  set  up  and  stuff  the  mesh. 

20feb87 

data  in  this  section  consists  of  groups  of  6  values  separated 
20feb87 

by  the  word  "cells’*,  those  values  whose  defaults  suffice 
20feb87 

need  not  be  input. 

20feb87 

20feb87 


input  word  definition 
20feb87 

20feb87 


(material  name)  (material  code)  material  name  abbreviation 
20feb87 

specifies  what  material  is  in  this  part  of 

20feb87 

the  mesh,  acceptable  names  are  air, tnt,btnt, 

20feb87 

pent , bpent , meth , bmeth , pbx , bpbx , anf o 


banfo, water.  20feb87 

=  code  must  agree  with  the  code  assigned  the 


20feb87 

20feb87 


material  in  the  sail  option  list 
(default  is  air=l) 


20feb87 

20feb87 

dx 

20feb87 

20feb87 

20feb87 

energy 

20feb87 

20feb87 


width  of  cells  in  this  part  of  the  mesh 
default  is  dx-0 . 0 

energy  density  to  be  inserted  in  cells  in 
this  part  of  the  mesh  (erg/gm) .  default 
is  ambient  energy  of  the  material. 


2Q£eb87 

20£eb87 

rho  mass  density  of  cells  in  this  part  of  mesh 
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20feb87 

20fab87 

20feb87 


(gm/cm**3) .  default  is  ambient  density. 


this 
1}  . 

the 


speed 

20feb87 

20feb87 

20feb87 

20feb87 

20feb87 

20feb87 

20feb87 

20feb87 

bound 

20feb87 


speed  of  material  in  this  part  of  the  mesh, 
at  least  the  first  cell  must  be  assigned  some 
small  initial  velocity  (if  less  than  umin, 
initial  velocity  will  not  be  present  by  cycle 
if  there  is  a  zero  velocity  in  the  first  cell 
program  will  stop, 
default  is  0.0  cm/sec. 

maximum  cell  number  for  which  the  previous 
five  values  hold  (default  -  bound=imax) 


20feb87 
20f eb87 
20£eb87 

these  six  words  may  appear  in  any  order  within  a  group. 
20feb87 

they  must  be  followed  by  the  word  "cells".  groups  must  be 


20feb87 

given  in  order  of  increasing  distance  from  the  center  of  the 


20feb87 


problem,  of  course  groups  may  not  be  mixed. 

20feb87 

20feb87 

the  final  word  of  data  input  on  cards  must  be  the 
20feb87 
word  "end". 

20feb87 

20feb87 

20feb87 

20feb87 

as  an  example,  a  data  set  for  a  new  problem  consisting 
20feb87 

of  100  cells,  75  pentolite  (the  first  cell  must  be  burned) 
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20febS7 

and  the  remaining  25  cells  air  at  ambient  conditions  follows. 
20feb8? 

fifty  stations  at  intervals  of  25.  cm  are  generated,  the 
20feb87 

material  codes  are  presumed  to  be  identical  with  those  input 
20feb87 

sail  options,  imax  and  nstat-no  of  stations  are  also  set  by 
20feb87 

sail  options. 

20feb87 

20feb87 

sap  prob=2001 . 00  t=l . e-4 

20feb87 

cycle-0. 0  input  cstop=2100.  detcel=1.0  rdw=11.0 
20feb87 

angle=0 . 0  head 
20feb87 


20feb87 
hob=0. 0 
2Qfeb87 
region 
20feb87 


sap  run  pentolite  100  cells 
stations  dx=25.  xs0=0 . 0 


bpent=2  speed=.01  dx=11.0  bound=l 
20feb87 


4  apr  75 


cells 


speed=.01  pent=3  dx=10.  bound=3  cells 
20feb87 

pent =3  bound=75  dx=10 . 0  cells 
20feb87 

air=l  dx— 3  0 .  bound=77  cells 

20feb87 

air=l  dx=42.  cells 
2QfebS7 


end 


20feb87 

20feb87 

20feb87 

20feb87 

deletion  of  the  statement  cycle=0.0  and  an  increment  of 
20feb87 

the  previous  run-s  stop  limit  (estop,  dstop,  or  ptstop)  with 


20feb87 


the  rest  of  the  data  unchanged  will  result  in  a  restart 
20feb87 

beginning  where  the  last  run  stopped,  any  changes  in  input 
20feb87 

will  take  precedence  over  tape4  values  obtained  from  the 
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20fefe>87 

=  the  zbiock,  except  of  course  those  parameters  which  change 

20feb87 

=  with  time  (eg  -  detcel,rdw, shkcel) .  if  restarting  before  a 

20feb87 

=  rezone  has  occurred,  the  inactive  portion  of  the  mesh 
20feb87 

=  (cells  shkcel+2  to  imax)  will  be  restuffed  with  the 

20feb87 

=  appropriate  data  from  the  "region"  section,  if  a  remesh  has 
20feb87 

=  taken  place  in  the  previous  run,  the  inactive  cells  are 
20feb87 

=  initialized  at  ambient  air  conditions  with  a  width  equal  to 
20feb87 

=  the  last  undisturbed  cell  in  the  mesh 

(x(shkceltl) -x(shkcel) ) .  20feb87 

~  in  this  case,  the  "region"  data  is  ignored  and  may  be  deleted 
20£eb87 

=  (omit  the  keyword  as  well) .  in  all  cases  "end”  must  mark  the 
20feb87 

*  end  of  the  sap  card  data. 

20feb87 

20feb87 

20feb87 

=  an  example  of  a  problem  using  sputter  data  as  input  is 

the  20feb87 
=  following. 

20feb87 

20feb87 

20feb87 

=  sap  t=i . e-3  prob=1313.00  cycle=0.00 

20feb87 

=  input  angle=0. 0 

20feb87 

-  head 
20feb87 

-  sputter  input  example$ 

20feb87 

=  sputin  yield=2.0  hob=5.e5  scale  fb=fb21  d»:~ .  2 

20feb87 
=  end 

20feb87 

20feb87 

=  where  the  sputter  tape  fb21  is  used  with  data  from  a  dump 

closest  2Qfeb87 

=  to  but  latter  than  l.e-3sec.  this  data  will  be  scaled  to  a 
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2  kt 


with 

card 


20feb87 

burst  at  5  km.  the  remaining  cells  will  be  .2  cm  filled 
20feb87 

ambient  air  conditions.  note  that  the  '$'  on  the  header 
20feb87 

is  for  disspla  plotting  routines  and  must  be  present 
20feb87 

as  a  string  terminator. 

20feb87 

20febS7 
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Appendix  B 


# 

#  Batch  File  jobsap.com 

0 

t  embedded  NQS  options 

0  @$-eo  0  directs  stderr  to  stdout 

0  @$-lt  0:30:00  0  per-process  CPU  time 
0  @$-lT  0:30:00  0  per-request  CPU  time 
0  Q$  0  end  of  embedded  NQS  options 

# 

0  Add  SAP  to  the  HULL  Library 
0 

cd  /u2/h5serag0/sap 

0 

In  /u2/h5serba0/HULL122/hulll22 .system  old 
0 

cat  >sail.in  «E0F 
sail  version  122.1  update 
*i  189999 
*b  runsap 

*keepto  runsap  runsap  eql 
♦read  modif.cr 
♦label  runsap 
EOF 

sail  <sail.in  >sail.s.out 
In  new  sapvl. system 
rm  -f  new  old 
0 

0  Re-Sequence  HULL  Library 
0 

In  sapvl. system  old 
0 

cat  >sail.in  «E0F 
sail  version  123  update 
seq  program  runsap 
EOF 

sail  <sail.in  >sail.s.out 
In  new  sapv2. system 
rm  -f  new  old 
0 

0  This  converts  the  program  runsap  from  the  hull  binary 
0  to  ascii  format 

0  (Stored  in  sapvl . system) 

0 

In  sapv2. system  old 
0 

cat  >sail.in  <<E0F 
sail  copy  convert 
EOF 
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sail  csail.in  >sail.s.out 
In  new  runsap. ascii 
rm  -f  new  old 
* 

i  This  prepares  coded  hull  library  of  runsap 

#  from  runsap. ascii 

# 

In  runsap. ascii  old 

* 

cat  >sail„in  «EOF 
sail  convert  copy 
EOF 

sail  <sail.in  >sail.s.out 
In  new  sapv3. system 
rm  -f  old  new 
# 

#  Now  Compile  and  Run  SAP 

# 

In  sapv3. system  old 

# 

cat  >sail.  in  «EQF 

sail  lineno  options  runsap=l  endoptions  program  sap 
EOF 

sail  <sail.in  >sail.s.out 
mv  source  sap.f 
# 

* 

cat  >segldr.in  «EOF 

preset=-indef 

abs=sap 

bin=sap.o 

EOF 

# 

cft77  -a  static  -1  sap.l  sap.f 
segldr  segldr . in 
# 

time  sap  <sap.in  >sap.out 
# 
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Appendix  F 


1 

sap  run 

prob 

pentolite 

bdry 

air opt 

1. OOOOE+OO 

bdry 

angle 

7 . 5000E+01 

bdry 

1*  0000E+02 

atmos 

2 

0. OOOOE+OO 

bsum 

0 

0. OOOOE+OO 

burnout 

0 

0. OOOOE+OO 

estop 

2 

0. OOOOE+OO 

CO 

1 

0. OOOOE+OO 

Cl 

5 

0. OOOOE+OO 

cycle 

0 

0. OOOOE+OO 

detcel 

.  0 
0. OOOOE+OO 

dstop 

1 

0. OOOOE+OO 

dt 

1 

0. OOOOE+OO 

dtmin 

1 

0. OOOOE+OO 

dxmul 

1 

0. OOOOE+OO 

eth 

8 

0. OOOOE+OO 

euler 

0 

0. OOOOE+OO 

geom 

3 

0.  OOOOE+OO 

grndr 

0 

0.  OOOOE+OO 

hfb 

0 

0. OOOOE+OO 

hob 

0 

100  cells 

z  block 

1. 1100E+02 
1. 000QE+00 
0 . 0000E+00 
. 0000E+00 
. 0000E+Q0 
. 0000E+00 
. 1000E+03 
. 8000E+00 
. 0000E"01 
.  0000E+00 
. OOOOE+OO 
. 0000E+09 
. 0000E-09 
. Q000E+20 
.0000E+00 
.  3418E+16 
.0000E+00 
.OOOOE+OO 
.0000E+00 
.OOOOE+OO 
.OOOOE+OO 


0. 0Q00E+00 

iiaax 

0 . G000E+0Q 

1 . 0000E+02 

Ibraf 

0 . 0000E+QC 

1.000GE400 

method 

0 . OQQOE+OO 

1 . 0000E+00 

nh 

0 . 0000E+00 

5. 0000E+00 

nstat 

0 . 0QG0E+00 

0. 0000E+00 

nvarpr 

0.0000E+00 

2 . 5000E+01 

nsblk 

0 . 0000E+00 

4 . 2000E+01 

prdv 

0.0000E+00 

0. 0000E+00 

ptstop 

0 . 0000E+Q0 

1 . 0000E+02 

qmin 

0 . 00Q0E+00 

1.0000E-01 

rad 

0 . 0Q00E+00 

0. 0000E+00 

radsum 

0 . 0000E+00 

0. 0000E+00 

rdw 

0 . 0000E+00 

1 . 1000E+01 

rezone 

0 . 0000E+00 

1. 0000E+00 

shkcel 

0 . 0000E+00 

4. 0000E+00 

stabf 

0.0000E+00 

8.0000E-01 

t 

0 . OOOOE+OQ 

1 . 0000E-04 

thm 

G.0000E+00 

7 . 0947E+08 

main 

0.0000E+00 

1 . 0000E-04 

vise 

0 . 00Q0E+00 

1. 0000E+00 

yield 

0 . 0000E+00 

0. 0000E+00 

nmat 

0 . Q000E+0Q 

3. 0000E+00 

bpent 

0 . OOOOE+OO 

1. 0000E+00 

pent 

0 . 0000E+00 

2. 0000E+00 

air 

0. 0000E+00 

3. 0000E+00 

sap  run 

pentolite 

100  ce 

29  july  92 
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2 


block: 


bpent 

1 

i 

1. iOGOE+Ol 

pent 

*•* 

75 

7 . S100E+02 

air 

3 

100 

1 . 7770E+03 

1  sap  run 

pentolite 

100  cells 

29  july  92 

problem  111.00000  cycle 
total  energy  eth  =  8.34180E+16 


0.  t  =  1 . 00000E-04 

total  mass  thm  =  7.09472E+08 


raesh  values  start  of  this  run 

j  3£  u  p  epg 

rho  temp  zm  dx  j 

1  1. 10000E+01  1 . OOOQOE-02  1.03766E+11  5.15500E+10 

1 . 66000E+00  0 . 00000E+00  2.20946E+03  1.X0000E+01  1 

2  2 . 100Q0E+01  1.00000E-02  1.01325E+06  1.00000E+08 

1 . 66000E+Q0  0 . 00000E+00  1.31638E+04  1.00000E+0X  2 

3  3.10000E+01  1.00000E-02  1.01325E+06  1.00000E+OS 

1 . 66000E+00  0 . 00000E+0G  3.40798E+04  1.00000E+Q1  3 

4  4.10000E+01  0 . OOOOOE+OO  1.01325E+06  1.00000E+08 

1 . 66000E+00  C. 00000E+00  6.49558E+04  1.0G0Q0E+01  4 

5  5 . 10000E+01  0. OOOOOE+OO  1.01325E+06  1.00000E+08 

1 . 66000E+00  0. OOOOOE+OO  1.05792E+05  1.00000E+C1  5 

6  6 . 10000E+01  0. OOOOOE+OO  1.01325E+06  1.00000E+08 

1 . 66000E+00  0. OOOOOE+OO  1.56588E+05  1.00000E+01  .  6 

7  7 . 10000E+01  0. OOOOOE+OO  1.01325E+06  1.00000E+08 

1 » 66000E+00  0. OOOOOE+OO  2.17344E+05  1.00000E+01  7 

8  8.10000E+01  0. OOOOOE+OO  1.01325E+06  1.00000E+08 

1 . 66000E+00  0. OOOOOE+OO  2.88060E+05  1.00000E+01  8 

9  9 . 10000E+01  0. OOOOOE+OO  1.01325E+06  1.00000E+08 

X . 66000E+00  0. OOOOOE+OO  3.68736E+05  1.00000E+01  9 

10  1 . 0100QE+02  0. OOOOOE+OO  1.01325E+06  1.00000E+08 

1 . 66000E+00  0. OOOOOE+OO  4.59372E+05  1.00000E+01  10 

11  1 . 11OO0E+O2  0. OOOOOE+OO  1.01325E+06  1.00000E+08 

1 . 66000E+00  0. OOOOOE+OO  5.59968E+05  1.00000E+01  11 

12  1 » 21000E+02  0, OOOOOE+OO  1.01325E+06  1.00000E+08 

1 . 66000E+00  0. OOOOOE+OO  6.70524E+05  1.00000E+01  12 

13  1 . 31000E+02  0. OOOOOE+OO  1.01325E+06  1.00000E+08 

1 . 66000E+00  0. OOOOOE+OO  7.91040E+05  1.00000E+01  13 

14  1 . 41000E+02  0. OOOOOE+OO  1.01325E+06  1.0Q000E+08 

1,660001+00  0. OOOOOE+OO  9.21516E+05  1.00000E+01  14 

15  1 . 51000E+02  0. OOOOOE+OO  1.01325E+06  1.00000E+08 

1 . 66000E+00  0. OOOOOE+OO  1.06195E+06  1.00000E+G1  15 

16  1 . 61000E+02  0. OOOOOE+OO  1.01325E+06  1.00000E+08 
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1.66000E+00  0 « OOOOOE+OO  1.2X235E+0S  l.OOOOOE+OX 

17  1.71000E+02  0. OOOOOE+OO  1.01325E+06  1 

1 - 66000E+00  0. OOOOOE+OO  1.37270E+06  l.OOOOOE+Oi 

18  1 » 81000E+02  0. OOOOOE+OO  1.01325E+Q6  1 

1 . 66000E+00  0, OOOOOE+OO  1.54302E+06  1.00000E+01 

19  1 . 91000E+02  0. OOOOOE+OO  1.01325E+06  1 

1 . 66000E+00  0 . OOOOOE+OO  1.72330E+06  1.00000E+01 

20  2 . 01000E+02  0. OOOOOE+OO  1.01325E+06  1 

1 . 66000E+00  0 «, OOOOOE+OO  1.91353E+06  1.00000E+01 

21  2 . 11000E+02  0. OOOOOE+OO  1.01325E+06  1 

1 . 660Q0E+00  0. OOOOOE+OO  2.11373E+06  1.00000E+01 

22  2 . 21000E+02  0, OOOOOE+OO  1.G1325E+06  1 

1 . 66000E+00  0. OOOOOE+OO  2.32388E+06  l.GOOOOE+OX 

23  2 • 31000E+02  0. OOOOOE+OO  1.01325E+06  1 

1 . 66000E+00  0. OOOOOE+OO  2.54400E+06  1.00000E+01 

24  2 . 41000E+02  0. OOOOOE+OO  1.01325E+06  1 

1 . 66000E+00  0, OOOOOE+OO  2.77408E+06  1.00000E+01 

25  2 . 51000E+02  0. OOOOOE+OO  1.01325E+06  1 

1 . 66000E+00  0. OOOOOE+OO  3.01411E+06  1.00000E+01 

26  2 . 61000E+02  0. OOOOOE+OO  1.01325E+06  1 

1 , 66Q00E+00  0. OOOOOE+OO  3.26411E+06  1.00000E+01 

27  2.7 1000E+02  0. OOOOOE+OO  1.01325E+06  1 

1 . 66000E+00  0. OOOOOE+OO  3.52406E+06  1.00000E+01 

28  2 . 81000E+02  0. OOOOOE+OO  1.01325E+06  1 

1 » 66000E+00  0. OOOOOE+OO  3.79398E+06  1.00000E+01 

29  2 . 91000E+02  0. OOOOOE+OO  1.Q1325E+06  1 

1 . 66000E+00  0. OOOOOE+OO  4.07386E+06  1.00000E+01 

30  3 . 01000E+02  0. OOOOOE+OO  1.01325E+06  1 

1 . 66000E+00  0. OOOOOE+OO  4.36369E+06  1.00000E+01 

31  3.11000E+02  0. OOOOOE+OO  1.01325E+06  1 

1 . 66000E+00  0. OOOOOE+OO  4.66349E+06  1.0Q0G0E+01 

32  3 . 21000E+02  0. OOOOOE+OO  1.01325E+06  1 

1 . 6600QE+00  0. OOOOOE+OO  4.97324E+06  1.00000E+01 

33  3 . 31000E+02  0. OOOOOE+OO  1.01325E+06  1 

1 . 66000E+00  0. OOOOOE+OO  5.29296E+06  1.00000E+01 

34  3 . 410OOE+O2  0. OOOOOE+OO  1.01325E+06  1 

1 . 66000E+00  0. OOOOOE+OO  5.62264E+06  1.00000E+01 

35  3 . 51000E+02  0. OOOOOE+OO  1.Q1325E+06  1 

1 . 6609GE+00  0. OOOOOE+OO  5.96227E+06  1.00000E+01 

36  3 . 61000E+02  0. OOOOOE+OO  1.01325E+06  1 

1 . 66000E+00  0. OOOOOE+OO  6.31187E+06  1.00000E+01 

37  3 . 71000E+02  0. OOOOOE+OO  1.01325E+06  1 

1.66000E+00  0. OOOOOE+OO  6.67142E+06  X.OOOOOE+Ol 

38  3 . 81000E+02  0. OOOOOE+OO  1.01325E+06  1 

1 . 66000E+00  0. OOOOOE+OO  7.04094E+06  1.00000E+01 

39  3 . 91000E+02  0. OOOOOE+OO  I.01325E+06  1 

1 • 66000E+00  0. OOOOOE+OO  7.42042E+06  1.00000E+01 

40  4 . 01000E+02  0. OOOOOE+OO  1.01325E+06  1 

1 . 66000E+00  0, OOOOOE+OO  7.80985E+06  1.Q0000E+01 

41  4 . 11000E+02  0. OOOOOE+OO  1.01325E+06  1 

1 . 66000E+00  0. OOOOOE+OO  8.20925E+06  1.00000E+01 

42  4 . 21000E+02  0. OOOOOE+OO  1.01325E+06  1 
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X.66000E+Q0  O.OOOOGE+OO  8.61860E+Q6  l.OOOGOE+Ol  42 

43  4 . 31000E+02  O.OOOOOE+OO  1.01325E+06  1.G0QO0E+08 

1 - 66000E+00  O.OOOOOE+OO  9.03792E+06  l.OOOOOE+Ol  43 

44  4 . 41000E+02  O.OOOOOE+OO  1.01325E+06  1.00000E+08 

1 . 66000E+00  O.OOOOOE+OO  9.46720E+06  l.OOOOOE+Ol  44 

45  4 . 5X000E+02  O.OOOOOE+OO  1.01325E+06  1.00000E+08 

X.66000E+00  O.OOOOOE+OO  9.90643E+06  X.OOOOOE+Ql  45 

46  4 . 61000E+02  O.OOOOOE+OO  1.01325E+Q6  1.00000E+08 

1.66000E+00  O.OOOOOE+OO  1.03556E+07  X.OOQOOE+Ol  46 

47  4 . 71000E+02  O.OOOOOE+OO  1.01325E+06  1.00000E+08 

X . 66000E+00  O.OOOOOE+OO  1.08148E+07  l.OOOOOE+Ol  47 

48  4 . 81000E+02  O.OOOOOE+OO  1.01325E+06  1.00000E+08 

X . 68000E+00  O.OOOOOE+OO  1.12839E+07  l.OOOOOE+Ol  48 

49  4 . 9X000E+02  O.OOOOOE+OO  1.01325E+06  1.00000E+08 

1 . 66000E+00  O.OOOOOE+OO  1.17630E+07  l.OOOOOE+Ol  49 

50  5 . 0X000E+02  O.OOOOOE+OO  1.01325E+06  1.00000E+08 

X . 66000E+00  O.OOOOOE+OO  1.22520E+07  l.OOOOOE+Ol  50 

X  sap  run  pentolite  100  cells  29  july  92 


problem  111.00000  cycle  *=  0.  t  =  1.00000E-04 

total  energy  eth  =  8.34180E+16  total  mass  thm  =  7.09472E+08 

mesh  values  start  of  this  run 

j  x  u  p  epg 

rho  temp  2m  dx  j 

51  5 . 11000E+02  O.OOOOOE+OO  1.01325E+06  1.00000E+08 

1 . 66000E+00  0 . 0000QE+00  1.27510E+07  l.OOOOOE+Ol  51 

52  5 . 21000E+02  O.OOOOOE+OO  1.01325E+06  1.00000E+08 

1 . 66000E+00  O.OOOOOE+OO  1.32600E+07  l.OOOOOE+Ol  52 

53  5. 31000E+02  O.OOOOOE+OO  1.01325E+06  1.00000E+08 

1 . 66000E+00  O.OOOOOE+OO  1.37789E+07  l.OOOOOE+Ol  53 

54  5 . 41000E+02  O.OOOOOE+OO  1.01325E+06  1.00000E+08 

1 . 66000E+00  O.OOOOOE+OO  1.43078E+07  l.OOOOOE+Ol  54 

55  5 . 51000E+02  O.OOOOOE+OO  1.01325E+06  1.00000E+08 

1 . 66000E+00  O.OOOOOE+OO  1.48466E+07  l.OOOOOE+Ol  55 

56  5 . 61000E+02  O.OOOOOE+OO  1.01325E+06  1.00000E+08 

1 . 66000E+00  O.OOOOOE+OO  1.53954E+07  l.OOOOOE+Ol  56 

57  5 . 71000E+02  O.OOOOOE+OO  1.01325E+06  1.00000E+08 

1 . 66000E+00  O.OOOOOE+OO  1.59541E+07  l.OOOOOE+Ol  57 

58  5. 81000E+02  O.OOOOOE+OO  1.01325E+06  1.00000E+08 

1 . 66000E+00  O.OOOOOE+OO  1.65229E+07  l.OOOOOE+Ol  58 

59  5. 91000E+02  O.OOOOOE+OO  1.01325E+06  1.00000E+08 

1 . 66000E+00  O.OOOOOE+OO  1.71015E+07  l.OOOOOE+Ol  59 

60  6 . 01000E+02  O.OOOOOE+OO  1.01325E+Q6  1.00000E+08 

1 . 66000E+00  O.OOOOOE+OO  1.76902E+07  l.OOOOOE+Ol  60 

61  6 . 11000E+02  O.OOOOOE+OO  1.01325E+Q6  1.00000E+08 

1. 66000E+00  O.OOOOOE+OO  1.82888E+07  l.OOOOOE+Ol  61 

62  6 . 21000E+02  O.OOOOOE+OO  1.01325E+06  1.00000E+08 

1 . 66000E+00  O.OOOOOE+OO  1.88973E+07  l.OOOOOE+Ol  62 

63  6 . 31000E+02  O.OOOOOE+OO  1.01325E+06  1.00000E+08 
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1 . 66000E+00  Q.OOQOOE+OO  1.95158E+07  X.OOOOOE+Ol 

64  6 . 41000E+02  O.OOOOOE+OO  1 . 01325E+06  1 

1.66000E+00  O.OOOOOE+OO  2.01443E+07  X.OOOOOE+Ol 

65  6 . 51000E+02  O.OOOOOE+OO  1.01325E+06  1 

1 . 66000E+00  O.OOOOOE+OO  2.07828E+07  l.OOOOOE+OX 

66  6.61000E+02  O.OOOOOE+OO  1.01325E+06  1 

1 . 66000E+00  O.OOOOOE+OO  2.14311E+07  X.OOOOOE+Ol 

67  6.7 1000E+02  O.OOOOOE+OO  1.01325E+06  1 

1 . 66000E+00  O.OOOOOE+OO  2.20895E+07  X.OOOOOE+Ol 

68  6.81000E+02  O.OOOOOE+OO  1.01325E+06  1 

1 . 66000E+00  O.OOOOOE+OO  2.27578E+07  X.OOOOOE+Ol 

69  6 . 91000E+02  O.OOOOOE+OO  1.01325E+06  1 

1 . 66000E+00  O.OOOOOE+OO  2.34361E+07  X.OOOOOE+Ol 

70  7 . 01000E+02  O.OOOOOE+OO  1.01325E+06  1 

1 . 66000E+00  O.OOOOOE+OO  2.41243E+07  1.00000E+01 

71  7 . 11000E+02  O.OOOOOE+OO  1.01325E+06  1 

1 . 66000E+00  O.OOOOOE+OO  2.48225E+07  X.OOOOOE+Ol 

72  7 . 21000E+02  O.OOOOOE+OO  1.01325E+06  1 

1 . 66000E+00  O.OOOOOE+OO  2.55307E+07  1.00000E+01 

73  7 . 31000E+02  O.OOOOOE+OO  1.01325E+06  1 

1 . 66000E+00  O.OOOOOE+OO  2.62488E+07  X.OOOOOE+Ol 

74  7 . 41000E+02  O.OOOOOE+OO  1.01325E+06  1 

1 . 66000E+00  O.OOOOOE+OO  2.69769E+07  1.00000E+01 

75  7 . 51000E+02  O.OOOOOE+OO  1.01325E+06  1 

1 . 66000E+00  O.OOOOOE+OO  2.77149E+07  1.00000E+01 

76  7 . 81000E+02  O.OOOOOE+OO  1.01325E+06  2 

1.22500E-03  2 . 87938E+02  6.46981E+04  3.00000E+01 

77  8 . 11000E+Q2  O.OOOOOE+OO  1.01325E+06  2 

1 . 22500E—03  2 . 87938E+02  6.98644E+04  3.00000E+01 

78  8 . 53000E+02  O.OOOOOE+OO  1.01325E+06  2 

1.22500E-03  2 . 87938E+02  1.06867E+05  4.20000E+01 

79  8 . 95000E+02  O.OOOOOE+OO  1.01325E+06  2 

1.225Q0E-03  2.87938E+02  1.17927E+05  4.20000E+01 

80  9 . 37000E+02  O.OOOOOE+OO  1.01325E+06  2 

1.22500E-03  2 . 87938E+02  1.29531E+05  4.20000E+0X 

81  9 . 79000E+02  O.OOOOOE+OO  1.01325E+06  2 

1. 22500E-03  2 . 87938E+02  1.41680E+05  4.20000E+01 

82  1 . 02100E+03  O.OOOOOE+OO  1.01325E+06  2 

1 . 22500E— 03  2 . 87938E+02  1.54373E+05  4.20000E+01 

83  1 . 06300E+03  O.OOOOOE+OO  1.01325E+06  2 

1.22500E-03  2 . 87938E+02  1.67610E+05  4.20000E+01 

84  1 . 10500E+03  O.OOOOOE+OO  1.01325E+06  2 

X . 22500E-03  2.87938E+02  1.81393E+05  4.20000E+01 

85  X . 14700E+03  O.OOOOOE+OO  1.01325E+06  2 

1.22500E-03  2 . 87938E+02  1.95719E+05  4.20000E+01 

86  1 . 18900E+03  O.OOOOOE+OO  1.01325E+06  2 

1 . 22500E— 03  2 . 87938E+02  2.10591E+05  4.20000E+01 

87  1. 23100E+03  O.OOOOOE+OO  1.01325E+06  2 

1.22500E-03  2 . 87938E+02  2.26007E+05  4.20000E+01 

88  1 . 27300E+03  O.OOOOOE+OO  1.01325E+06  2 

1 . 22500E— 03  2 . 87938E+02  2.41967E+05  4.20000E+01 

89  1 . 31500E+03  O.OOOOOE+OO  1.01325E+06  2 
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1.22500E-03  2 . 87938E+02  2.58472E+05  4.20000E+01 

90  1.35700E+03  0.000Q0E+00  1.01325E+06  2 

1 , 22500E-03  2.87938E+02  2.75521E+05  4.200Q0E+01 

91  1 . 39900E+03  0.00000E+00  1.01325E+06  2 

1 , 22500E-03  2.87938E+02  2.93115E+05  4.20000E+01 

92  1 . 44 1G0E+03  O.OOOOOE+OO  1.01325E+06  2 

1 . 22500E-03  2.87938E+02  3.11254E+05  4.2000QE+01 

93  1.48300E+03  O.OOOOOE+OO  1.01325E+06  2 

1.22500E-03  2 . 87938E+02  3.29937E+05  4.20000E+01 

94  1 . 52500E+C3  O.OOOOOE+OO  1.01325E+Q6  2 

1.22500E-03  2 . 87938E+02  3.49165E+05  4.20000E+01 

95  1 . 56700E+03  O.OOOOOE+OO  1.01325E+06  2 

1.2250OE-O3  2 . 87938E+02  3.68937E+05  4.20000E+01 

96  .  1 . 60900E+03  O.OOOOOE+OO  1.01325E+G6  2 

1 . 22500E-03  2 . 87938E+02  3.89254E+05  4.20000E+01 

97  1 . 65100E+03  O.OOOOOE+OO  1.01325E+06  2 

1 . 22500E-03  2 . 87938E+02  4.10115E+05  4.20000E+01 

98  1 . 69300E+03  O.OOOOOE+OO  1.01325E+06  2 

1.22500E-03  2 . 87938E+02  4.31521E+05  4.20000E+01 

99  1 . 73500E+03  O.OOOOOE+OO  1.01325E+06  2 

1 . 22500E-03  2 . 87938E+02  4.53472E+05  4.20000E+01 

100  1.77700E+03  O.OOOOOE+OO  1.01325E+06  2 

1.22500E-03  2 . 87938E+02  4.75966E+05  4.20000E+01 

1 trouble  in  subroutine  burn  -  abort 
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